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Abstract 

The two-photon Rabi Hamiltonian is a simple model describing the interaction of light with matter, 
with the interaction being mediated by the exchange of two photons. Although this model is exactly 
soluble in the rotating-wave approximation, we work with the full Hamiltonian, maintaining the 
non-integrability of the model. We demonstrate that, despite this non-integrability, there exist 
isolated, exact solutions for this model analogous to the so-called Juddian solutions found for the 
single-photon Rabi Hamiltonian. In so doing we use a Bogoliubov transformation of the field 
mode, as described by the present authors in an earlier publication. 

PACS number(s): 03.65.-w, 42.50.-p, 32.80.-t 



I Introduction 



The Rabi Hamiltonian (RH), introduced by Rabi in 1937 |jj , has long served as a popular and successful 
model of the interaction between matter and electromagnetic radiation. The Hamiltonian provides 
a description of an atom approximated by a two-level system interacting via a dipole interaction 
with a single mode of radiation . Typically, this Hamiltonian is studied within the rotating- wave 
approximation (RWA), which results in the well-known Jaynes-Cummings model (JCM) j|. The JCM 
is exactly integrable, whereas the full RH is not. 

The two-photon Rabi Hamiltonian (TPRH) is an obvious extension of the original RH, where the 
atomic transitions arc induced by the absorption and emission of two photons rather than one. Such 
non- linear optical processes have been of considerable interest [4], with applications including two- 
photon lasers and two-photon optical bistability [5]. The TPRH is not known to be integrable, whereas 
its RWA counterpart is, as has been demonstrated by Sukumar and Buck [6] and Compagno and Persico 
[7]. It should be noted from the outset that the TPRH, and its RWA variant, are phenomenological 
Hamiltonians, in that they neglect the effects of intensity-dependent Stark shifts of the atomic levels ||). 
Nevertheless they do provide useful prototypes of two-photon interactions |5), and their similarity 
with the RH and its RWA variant allows fruitful comparisons to be made fllPf , The TPRH has 
been of considerable theoretical interest due to the connection of the two-photon interaction to the 
group SU(1, 1) and to the squeezed states ]H[| . Experimentally, the observation of two-photon Rabi 
oscillations in experiments with Rydberg atoms [jl2| , |l3|] has also contributed to the interest in this type 
of model. 

Comparatively little attention has focused on the TPRH without the RWA. A notable exception to this 
is the work of Ng et al. |l4| , who have used numerical diagonalisation in a truncated basis to investigate 
the spectrum and simple dynamics of the system. Their analysis indicates that the spectrum of the 
full Hamiltonian is significantly different to the RWA spectrum, and that making the RWA also alters 
appreciably the dynamics. These results fit in well with other work regarding the RWA, which suggest 
that the consequences of making this approximation may be greater than usually thought [jl5],[l(| . Here 
we shall exclusively consider the full Hamiltonian without the RWA. 

In this communication we discuss a number of exact results for the TPRH without the RWA. After 
introducing the Hamiltonian and examining its connection to the group SU(1, 1), we consider the limit 
of the Hamiltonian in which the atomic levels become degenerate. In so doing we obtain a condition on 
the range of atom-field couplings for which this model remains mathematically valid. We then proceed 
to obtain a set of isolated, exact solutions for the Hamiltonian. Their counterparts are well known for 
the single-photon RH, where they are referred to as Juddian solutions [jl7|,[l8). Such solutions tell us a 
great deal about the structure and symmetries of this type of non-adiabatic model. They may also serve 
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as benchmarks for numerical techniques and as foundations for perturbative treatments. Furthermore, 
the existence of isolated exact solutions in non-integrable quantum models is also of interest from the 
perspective of studying possible quantum chaos in such systems [jl^,^o|. In determining these solutions 
for the TPRH we shall utilise an appropriate Bogoliubov transformation of the field mode, an approach 
outlined by the present authors in a previous publication pl|, to be referred to as I hereafter. 



II The Hamiltonian 



The TPRH describes the interaction of a two-level atom with a single bosonic field mode via a two- 
photon interaction. The field is described by the annihilation and creation operators, b and re- 
spectively, which obey the usual commutation relation, [b,b'~\ = 1. The two-level atom is described 
by the Pauli pseudo-spin operators {c^fc = x,y,z}, which satisfy the SU(2) commutation relations, 
[cr x , <j y ] = 2ia Zl plus cyclic permutations. We define the raising and lowering operators as a+ = (Tx+iay, 
<r_ = <t~ — ia„. 



In terms of these operators, the TPRH is given by 

H 2 , = y a, + ootfb + g ((6+) 2 + b 2 ) (a + + , 



(1) 



where loq is the atomic level splitting, to is the frequency of the boson mode and g is the coupling 
strength of the atom to the field. Note that here we have the operators (b^) and b 2 inducing atomic 
transitions, instead of W and b, as we would have in the single-photon RH. It is convenient to rescale 



the Hamiltonian as H- 



2 7 



LjH2-f, where 



H 2l = iua z 



tfb + \((tf) 2 + b 2 )a x , 



(2) 



and Co = ^ 



and A = ^S.. The TPRH is not known to be integrable. Like the single-photon RH, the 
TPRH possesses a conserved quantum number, which in the present case is 7T2 7 , namely an eigenvalue 
of the operator 
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exp 
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2 

cr z exp 



(b+6- 



1 



-6+6 



(3) 
(4) 



It is simple to show that [i?2 7 , n2 7 ] = 0. The operator exp (i^tfb) is the square root of the elementary 
parity operator exp (iirtf&J = H 2 ^, and has been denoted as the Fourier operator. Its role in squeezing 
has been described in detail elsewhere p2| . 

As the level- flip in the TPRH is induced by two photons, the condition for the system to be on 
resonance is loq — 2w, or alternatively that Co = 1. 
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II.l The TPRH and the 377(1,1) Group 

The TPRH contains only quadratic combinations of the bosonic creation and annihilation operators. 
Consequently, we may write the Hamiltonian in terms of the three operators Ko, K + and K-, which 
are defined as 

k^ 1 -^) 2 , K^K\ = \b\ jsr = I 6 t 6+ i. (5) 

These operators form a closed Lie algebra SU (1, 1), defined by the commutator relations 

[K ,K ± }=±K ± , [K.,K+]=2K . (6) 

The corresponding invariant Casimir operator C is given by 

C^Kl-±(K + K_+K_K + ), (7) 

which commutes with all three generators. For our purposes here, we shall use a unitary irreducible 
representation of this algebra known as the positive discrete series D + (fc) In this representation 
the basis states {\k,m}} diagonalise the operator Kq 

K \k,m) = (m + k)\k,m), (8) 

for k > and m = 0, 1, 2, . . .. The action of the Casimir operator in this representation is 

C\k,m) = k(k-l)\k,m). (9) 

The operators K + and K- are Hermitian conjugate to each other and act as raising and lowering 
operators respectively within D + (k), 



K + \k,m) = \J(m + 1) (m + 2k)\k, m + 1), 



K-\k,m) = (m + 2k- l)\k,m- 1). (10) 

For the single-mode bosonic realisation of SU(1, 1) that we require here, the Bargmann index k is 
equal to either j or |. In terms of the original Bose operators the states \k, m) are given equivalently 
as 

|1/4 ' m> "(i)!^ 2 " l|0); |3/4 ' m> "(2^TTj! (6t)2m+1|0> ' m = ' 1 ' 2 '"- ^ 
So by switching to a SU (1,1) representation we explicitly acknowledge that we are splitting the Hilbert 
space of the boson field into two independent subspaces. Each subspace is labeled by the Bargmann 
index k — j, | and only contains either all even (k — j) or all odd (k = |) number states . It is 
interesting to note in passing that the algebra SU(l, 1) may also be used to describe a system of two 
bosonic modes, which interact in such a way as to preserve the total particle number [^5|^4| 
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In terms of the SU(1, 1) generators, the TPRH can be written 



H 2 ^ = t ^a z + 2cj(K -~)+2g(K + +K_){a + +a-), (12) 



2 * V 4 
with the rescaled Hamiltonian being given by 



H 2l = uht z + 2 (k - i \ + 2A (K+ + K-) o x . (13) 
II.2 Squeezing and 377(1,1) 

The relationship between the group SU(1, 1) and squeezing has been described in detail elsewhere |2(| 
and the use of squeezed states in finding exact isolated solutions has been previously discussed in I. 
Here we simply note that the general squeezing operator S can be written as 

S (a, (3) = exp (aK+) (l - | CT | 2 )^ exp (-a*K-) exp ((3 (K - 1/2)) , (14) 

where a and (3 are squeezing parameters, with (3 real and a a complex number with modulus \a\ < 1. 
S is a unitary operator, S'S = 1, and provides a representation of the group SU(1, 1). With it we 
may make unitary transformations of the bosonic annihilation and creation operators, such that 

S{a,/3)bSHa,/3) = e~ 10 (l - \<r\ 2 y 1/2 (b - atf) = c, 

S{? t p)tf&{o,p) = e 10 (l-\a\ 2 y 1/2 (tf -a*b) =c f . (15) 

The operators c and c' satisfy the commutation relation [c, <7] =1, and are henceforth referred to as 
squeezed boson operators. 

Ill Degenerate atomic levels 

For degenerate atomic levels, ujo = = ui, the (rescaled) TPRH takes the form 

H(;°= 0) =tfb + \((tf) 2 + b 2 )a x . (16) 

Consequently eigenstates of H^ 0-0 " 1 are also eigenstates of a x , and we are led to consider the bosonic 
Hamiltonian, 

/47= 0) =&t&±A((&t) 2 + 6 2 ), (17) 

where the two signs correspond to the two eigenvalues of <r x . The Hamiltonian of Eq. ( |T7| ) has the form 
of a squeezed harmonic oscillator. In seeking its eigen-solutions, it is convenient to use the squeezed 
bosons discussed above. 
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Inverting the relations (15), setting [3 to zero and constraining a to be real, we obtain the following 
forms for the squeezed bosonic operators 

t tf±ab b±atf 

ci s 77=T>- c± s 7T=7i- (18) 



where the subscript on these operators corresponds to the sign in Eq. (p_7[). We now choose a to be 
given by 

<r = ^; fi = VT^. (19) 

Writing the Hamiltonian in terms of these squeezed operators with this value of a we have 

h^= 0) = « ( ct ±C± + 1) (20) 

The eigenstates of this Hamiltonian are clearly the number states of the c±-type bosons, which we 
denote |n, =f<t), such that c±c±\n,^fa) — n\n,^f<r). In our original unsqueezed representation, these 
states have the form 



\n;T<r) = — 



< 1/4 'tf±ab 



e^ bt2 |0), (21) 



showing them to be the usual squeezed number states B6[ , Thus we see the eigenenergies of the 
Hamiltonian of Eq. ( po|) to be 

E(-o=o) = | ri+ lj n _l ) (22) 
such that h<£ 0=0) \n; Tct) = ^ 0= °V; To)- 

An important feature of the TPRH is revealed by considering this ujq = case. As we saw above, the 
eigenfunctions of the bosonic part of the loq = Hamiltonian are number states of the squeezed bosons. 
The squeezed vacuum, |0; T^ 7 ), is proportional to exp f T i CT &t 2 ^ | ) and this state is only normalisable 
for \a\ < 1, which corresponds to the conditions 

\M<h I77I <!• (23) 

Above this value of g/ui the uiq — Hamiltonian does not possess normalisable eigenfunctions and is 
thus unphysical. As has been discussed by Ng et al. fli}| , and as is borne out by numerical diagonali- 
sation |27|| , this analysis still holds for the loq 7^ case as the remaining operator in the Hamiltonian, 
uj(T z , is clearly a bounded operator (^8|, and thus presents no problems. Thus we see that the TPRH 
is only well defined for values of |A| < 1/2. This restriction on the coupling is not a severe restriction 
as, at higher couplings, effects not included in this model, such as the contributions of off-resonant, 
one-photon processes, will come into play, thereby compromising the physical relevance of the model. 
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IV Isolated exact solutions 



We now demonstrate the existence of a class of isolated, exact solutions for the TPRH, similar to the 
Juddian solutions found for the one-photon RH. Following I we seek solutions by first performing a 
Bogoliubov transformation of the field mode. Bearing in mind the to = result, we shall make the 
transformation from the original bosons b and w to the squeezed bosons c and c' , 



1/2 



(24) 



where we have assumed that (3 — and that a is real and to be determined. The justification of this 
is provided by subsequent results. 

With this change in bosonic representation, the rescaled TPRH becomes 

1 



H 2l = OJG z 



1 \* 

— I ere + crc 1 

K 

itt 1 + ■'>(■ 



1 + cr 2 ) c)c + a' 



4<7C t C+ 2(7) G x , 



(25) 



where k = (l — er 2 ). We now use an appropriate matrix representation for the Pauli matrices, which 
for this model is one in which a x is diagonal. We use 



(26) 



In terms of the two-component wavefunction, l^) = (j^j), the time-independent Schrodinger equation 
for the system, H2 1 \^) = E\^}, where E = E/lo is the rescaled energy, then reads 





1 







i 
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a x = 





















-1 


-i 




1 






cD|*i) 



\{[a + \(l + a 2 )] (c^ + c 2 ) 

+ [1 + a 2 + 4Aa] c f c + a 2 + 2\a - ke} = 0, 

i{[a-A(l + . 2 )] (c^ + c 2 ) 

+ [l + cr 2 - 4Act] c f c + a 2 - 2Ao- - kE\ |* 2 ) = 



(27) 



(28) 



It is immediately clear that if we set either [cr + A (l + cr 2 )] or [cr — A (l + cr 2 )] equal to zero, we make 
a determination of a and reduce either Eq. (|27]) or Eq. (^) considerably. Choosing the first of these 
options, we have 

cr + A(l + cr 2 ) =0 (29) 



which gives 



■i±n 



(30) 
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where CI is as defined previously in Eq. (19). In order that \a\ < 1, we must henceforth choose the 
positive sign in Eq. (|30|), so that a — ► as A —> 0. Thus our squeezing parameter is 

Q - 1 



2A 



(31) 



Note that had we pursued the other option and set [cr — A (l + er 2 )] = 0, we would have obtained the 
same determination of the squeezing parameter as for the lu = case given by Eq. dl9| ) . This second 
solution will be discussed later. Proceeding with this choice of squeezing given by Eq. (|3l|), Eqs. ( p7|) 
and (E8h become 



ne'e- 



E 



1 Cl 

2~ 2 



|*i)=0, 



Jl 2 (V 2 +c 2 ) + (2 -SI 2 ) 



2^1 „t 



c 1 c 



+- (i - n) (2 + fi) - Q£ } |* 2 ) = o. 



For j^i) and l^) we now choose simple Ansatze in terms of the squeezed number states; 

N 

71=0 

N-2 

|*2> = ^2 q m |m;er), 



m=0 



which gives us the equations 

N-2 

o)g m |m;cr) 

m=0 

AT 

^ wp„|n;er) 

n=0 



n=0 
jV-2 







SI" 




jnQ- 




2" 


} 



|n; cr) = 0, 



5m ^ {-V^ - ^ 2 [vV + 1) (m + 2)|m + 2; a) 



+Vm (rn - l)|m - 2; cr) + (2 - SI 2 ) m|m; cr) 
' f (l - Cl) (2 + Cl) - 2CIE\ \m; a) \ = 



For the first equation to hold, we must have 







Cl~ 




|aS1- 


M- 


2 


} 



= 0. 



As pn 7^ by Ansatz, we must have SI A — (^E + i 
baselines as 

1 



(32) 
(33) 

(34) 
(35) 



(36) 



(37) 



(38) 



£2 

2 2 



E = -- + [N + -\ n 



0. This establishes the so-called energy 

(39) 
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Equating the remaining coefficients of the number states in Eqs. (36) and ( |37| ) gives us the following 
set of equations 



u>q m 

OJPn 



+ (to — N) ftp 

Vi -ft 2 



0; 



< to < N 



ft 



y/n(n- l)q n - 2 - 



Vi -ft 2 
o 



y/(n+l)(n + 2)q n 



ft"' 7 ' 



[n + 1 - (2n + N + 1) ft 2 ] g„ = 0; 



^2 



< n < N. 



(40) 



(41) 



where q n = for n < and rt > N — 2. From the second of these equations we see that the Hamiltonian 
in this squeezed representation only couples number states that are different by multiples of two (for 
example it couples \n; a) to \n + 2; a) and \n + 4; a) but not to |n + 3; a)). Therefore our Ansatze need 
only include either all even number states or all odd number states. This is equivalent to restricting 
the solutions to a sector of the Hilbert space with a given Bargmann index, k. We also see that the 
minimum value of N is 2. When N is even, we obtain N + 1 equations for the N + 1 coefficients 



Pa 



0,2,. 



0, 2, . . .N — 2. When N is odd, we obtain N equations for the N 



coefficients p n ; n = 1 , 3, . . . , N, q m ;m = 1, 3, ... N — 2. Requiring that the determinants of these 
equation sets are equal to zero gives us the compatibility conditions which establish the location of the 
Juddian points on the energy baselines. These conditions for the first three values of N are: 



2 -6ft 2 
6 - 1 Oft 2 







-up- = 



8 (3 - 30ft 2 + 35ft 4 ) + 2 (7 - f 7ft 



2 ^ 2 



-cD 4 



= 



N = 2, 
N = 3, 
N = 4. 



(42) 



For N even we have a polynomial of the order ^ in w 2 and ft 2 , and hence A 



corresponding polynomial is of order ( N 2 1 ) . The equations are symmetric under w 



For N odd, the 
—u) or A — > —A, 

as expected. The solutions of these polynomials locate the exact solution in cD— A space. 
As in the one-photon case, there exists another degenerate solution at every Juddian point. This 
solution can be found by placing \a — A (f + cr 2 )] = rather than \a + A (l + cr 2 )] = into Eqs. ( pjj ) 
and (|28|). This simply interchanges the equations and hence the roles of \&i and ^2- 



V Results and Discussion 

Solving the above complementary conditions we have calculated twelve Juddian points for the resonant 
TPRH with 2uj — u> = 1, corresponding to values of of 2 < < 7. These are displayed in Table [j], 
listed to the first fO significant figures, and we have used the original units of Eq. (|f|). The N = 2 
and = 3 cases have such simple complementary conditions that closed analytic expressions may be 
found. On resonance these are given by g = — 4= and E = i("7f — l) f° r N — 2, and g — | J and 
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E = j y7y jq — lj for N = 3. This set of twelve Juddian points is indicated on the energy schema of 
the Hamiltonian in Fig. [l], where the schema was obtained by approximate numerical diagonalisation 
via a standard configuration-interaction method, using a basis size of the lowest 501 harmonic oscillator 
states p7j . The energy baselines, E = \ {~\ + (iV + |) are also plotted. 

As is clearly seen from Fig. |l|, the Juddian solutions found by this method occur at level-crossings 
in the energy schema, but that they do not cover every crossing, as is the case with the one-photon 
Hamiltonian. Considering the quantum numbers 7T2 7 of the two intersecting lines at each crossing, we 
see that the above type of solution can only describe the crossings of states having 7T2 7 = +1 with 
ones having 7T2 7 = — 1, and of crossings of states having 7T2 7 = +i with ones having 7T2 7 = — i. The 
remaining four types of possible crossings are not described. This situation is summarised in Table ^. 
This series of crossings can be understood by considering the parity operator Hf 

n| 7 = exp (iirtfb) , (43) 

which obviously commutes with the Hamiltonian. From considering the eigenvalues of this operator 
we see that the Juddian solutions we have described occur between levels which have the same value of 
n§_. Thus although the Ansatze of the Juddian solutions above are not eigenstates of the Fourier-like 
operator n2 7 , they are eigenstates of parity. 

These crossings can be viewed in another way. Ng et al. have introduced a unitary transformation 
which decouples the spin and bosonic degrees of freedom [jl4| . After the application of this transfor- 
mation, the bosonic part of the Hamiltonian is given by 

h (M) = l Mwo (_i)*o-* + 2uj (k q --^j+Ag (K + + K_) , (44) 

where the numbers M = ±1, corresponding to the spin degree of freedom, and k — 4, |, corresponding 
to the Bargmann index, serve to characterize the four independent sub-spaces into which the full 
Hilbert space of the Hamiltonian decomposes under this transformation. We thus see that the crossings 
detailed above occur between states with the same Bargmann k index, but with different M indices. 
By contrast the missing crossings occur between states with different values of k. 
The reason why the above Ansatze can describe these solutions and not the others is as follows. The 
solutions that we have been able to find occur at crossings between energy eigenf unctions that both 
have the same Bargmann index k, which means that both states are composed of either all even or all 
odd number states. At the Juddian points these two eigenstates become degenerate in energy and thus, 
to find the energy at the level-crossing we may form a linear superposition of the two eigenstates, which 
will, in general, not be an eigen-state of n| 7 - Because the degenerate energy eigenstates belong to 
the same fc-sector, the formation of the superposition allows the individual terms in one wavefunction 
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to add to the terms in the other. If we form the superposition correctly, the resultant wavefunction 
may have a form much simpler than the constituent wavefunctions. This is exactly the case when we 
choose the Ansatz (|35l). 

The solutions that we have been unable to find with the above method occur at the level-crossings 
between energy eigenstates that have different Bargmann indices. This means that one eigenstate is 
composed of odd number states, whilst the other is composed only of even number states. Consequently, 
no superposition of these states will lead to a reduction in the complexity of either wavefunction and 
we have been unable to find simple Ansatze at these level-crossings. 

Although the above method is not directly extensible to the remaining crossings, it may still be the 
case that exact solutions can be found. Although there is no a priori reason to expect that these exact 
solutions exist, by looking at the energy schema generated numerically, we observe that the remaining 
level-crossings appear to lie along base-lines described by 

E = ~+nQ, n = 2,3,..., (45) 

where, as previously, SI = \/l — 4A 2 . These baselines are so similar to the baselines for the Juddian 
solutions found above that it would seem likely that Juddian solutions could also be found at these 
remaining crossings. 

VI Conclusions 

We have shown that a set of isolated, exact solutions exists for the two-photon Rabi Hamiltonian. 
These are seen to occur at a sub-set of the level-crossings in the energy schema and although we 
have not described every level crossing, we have been able to explain this in terms of the symmetry 
properties of the crossing states. 

This type of isolated solution also occurs in the single-photon Rabi Hamiltonian, where they are referred 
to as Juddian solutions. Several methods have been proposed for finding these solutions ^9|. In 
I we have given a method for finding these solutions that is directly analogous to the one used here, 
except that for the one-photon case we have used displaced, rather than squeezed, bosons. In this 
case, the Ansatz does provide solutions at every level-crossing in the spectrum. It should be noted 
that the single-photon Hamiltonian is a simpler model than the TPRH, as the conserved quantum 
number analogous to 7T2 7 only has two eigenvalues, and thus there is only one type of level crossing, 
whereas in the current model has four eigenvalues and there are six different types of crossing. 
It is hoped that the results presented here, in conjunction with those in I, will be of use in the analysis 
of this kind of non-adiabatic model. These isolated solutions seem an ideal starting point for the 
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analysis of this kind of situation, since they provide exact results which may be used as bench-marks 
for further methods. They also provide crucial insight into the symmetries of the model. 
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g 

Figure 1: The first twelve Juddian points (diamonds) of the two-photon Rabi Hamiltonian determined 
by the method outlined in the text, plotted against the energy spectrum determined numerically (solid 
lines). Also plotted are the energy baselines (dotted lines). Each point is labeled with its order N. 
The Hamiltonian is resonant; 2u> = cJq = !• 
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Table 1: The couplings (g), energies (E), as well as order (N)oi the first twelve Juddian points of the 
resonant two-photon Rabi Hamiltonian (2w = coq = !)• 
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Table 2: Description of the parities of the level-crossings described by the TPRH Juddian solutions 
presented here, "y" denotes that the crossing is described, "n" that it is not, whereas "-" indicates 
that no such crossing occurs. 
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